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Abstract. Wo study the random Fibonacci sequences defined hy Fi = F2 = Fi = F2 = 1 
and for n > 1, F„^2 = Fn+i i F„ (linear case) and F„^2 = l^n+i i Fn\ (non-linear case), 
where each ± sign is independent and either + with probability p or — with probability 1 — p 
(0 < p < 1). Our main result is that the exponential growth of Fn for < p < 1, and of Fn for 
l/3<p<lis almost surely given by 



JO 



log X duot{x), 

'0 

where a is an explicit function of p depending on the case we consider, and i/a is an explicit 
probability distribution on R_|_ defined inductively on Stern-Brocot intervals. 

In the non-linear case, the largest Lyapunov exponent is not an analytic function of p, since 
we prove that it is equal to zero for < p < 1/3. We also give some results about the variations 
of the largest Lyapunov exponent, and provide a formula for its derivative. 



1. Introduction 

In this article, we wish to investigate the exponential growth of random Fibonacci sequences 
{Fn)n>i and (F„)„>i, defined inductively by Fi = F2 = Fi — F2 — 1, and for all n > 1, 

(1) Fn+2 = F„+i ± F„ (linear case), 

(2) Fn+2 = \Fn+i ± K\ (non-linear case), 

where each ± sign is independent and either + with probability p or — with probability 1 — p 
(0 < p < 1). In the case p = 1/2, (|fra|) and (F„) have the same distribution law as the 
sequence (|in|) studied by Viswanath [TU]. In his paper, using Furstenberg's formula [1] (see also 
[1], Chapter II), Viswanath proves that with probability 1, 

1~\ > 1.13198824..., 

n — ' + 00 

where the logarithm of the number is computed as the integral of the function 

1 , / 1 + 4m' 
■log 



4 ""^ V(l + '^^)^ 

with respect to some explicit "fractal" measure vj. 

Our purpose here is to give a formula for any parameter p e]0, 1], and to provide some results 
on the dependence on p of the upper Lyapunov exponents. By contrast with the linear case ([T]), 
the non-linear case Q cannot be viewed as a product of i.i.d. random matrices. This explains 
why the upper Lyapunov exponent in the non-linear case is not an analytic function of p, as can 
be seen in Theorem 11.11 Our method does not make use of Furstenberg's formula, but relies on 
the reduction of random Fibonacci sequences exposed in [8] . 
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1.1. Results. Our main result is the following. 

Theorem 1.1. 

• Linear case: 



1 

-log|F„| > 7p = / logx di/a{x) a.s., 

n n^oo 



where 



3p-2 + J5p^ - 8p + 4 
a := 



2p 

• Non-linear case: If p < 1/3, — logF„ > 7p = 0. Ifp> 1/3, 



n 



n 

where 



1 ~ ^ f°° 

-logF„ > 7p = / \ogxdva{x) a.s., 

n n-*oo 



2p 



p + ^p(4 - 3p) 

In boti cases, Va is an explicit probability distribution on R_|_ dehned inductively on Stern-Brocot 
intervals (see Section \3~S and FigurelTji. 
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Figure 1. The measure Va on Stern-Brocot intervals of rank 1, 2, 3. First assign 
mass 1 — a to [0,1] and a to [l,oo]. Once i^a is defined on some Stern-Brocot 
interval of rank r, a proportion a of its mass is given to the left (respectively 
right) subinterval of rank ?' + 1 when r is odd (respectively even). 

For p — 1/2, we get a — (f>^^ both in the linear and the non-linear cases, where (f) :— (1 + \/5) /2 
is the golden ratio. The measure is then equal to Viswanath's fractal measure conditioned on 
R"*". For p = 1, which corresponds to the classical Fibonacci sequence, a = 1 and j/q. is the Dirac 
mass on 0. When p in the linear case, or p — > 1/3 in the non-linear case, a — > 1/2 and 
Va which is the probability measure on giving the same mass 2~'' to each Stern-Brocot 

interval of rank r. This measure is related to Minkowski's Question Mark Function (see [3]): 

Va;e[0,l], l{x) = 2y^,2{[Q,x]). 

Remark 1.2. The exponents jp and 7p correspond to almost-sure exponential growth of random 
Fibonacci sequences. We could also consider the average point of view, that is, ask for the limit 

of — log(E(|F„|)) and — log(E(F„)) (where E stands for the expectation), 
n n 

In the non-linear case, we know how to give an explicit expression of the limit for any p. (Of 
course, by Jensen's inequality, this limit is bounded below by^p.) It turns out that the critical 
value ofp for which this limit vanishes is p = 1/4 (compare it with the value 1/3 in Theorem ll.lfl . 



HOW DO RANDOM FIBONACCI SEQUENCES GROW? 



3 



The techniques used to obtam this result are quite different from those presented here; They 
mainly rely on ideas introduced in [5], in which the case p = 1/2 is treated. Details and proofs 
will be given in a forthcoming paper. 

In Section 14.31 wc study some properties of the functions p i— ^ 7p and p i— *■ 7p, and prove the 
following theorem. 

Theorem 1.3. 

• Linear case: p ^ is an increasing function of p, satisfying 

l™7p = 0, 7i=log(/) and = i^^!!^. 

p-»o dp 2 

• Non-linear case: p >-^jp is a continuous function of p, increasing on ]l/3, 1], satisfying 

> = 0forp<l/3, 7i=log0 and ^(i) ^ i2|^. 

dp 2 

One of the ingredients for the proof is a formula for the derivative of 7p (or 7^) with respect to 
a, involving the product measure Va ® (see Proposition 14. 5p . 



7p 

f > 7(a) A 




1/3 1/2 I p 1/2 1 a 



Figure 2. Graphs of 7p and 7^ as a function of p (left) and as a function of a 
(right). Notice that 7p = for p e]0, 1/3]. For p e [1/3, 1], the graph of 7^ is not 
the straight line it looks like: Compare the average slope with the derivative in 1. 



2. Reduced sequences 

2.1. Random paths in T and T. The sequences {Fn)n and {Fn)n can be read along random 
paths in the trees T and T described as follows. These two trees have the same structure, but 
differ by the labels attached to the vertices. Each vertex is labelled by an integer: The root and 
its only child are labelled by 1. Any other vertex has two children, left and right. If tj is a right 
child, its label is the sum of its father's and grandfather's labels; If w is a left child, its label is 
the difference between its father's and grandfather's labels in the tree T, whereas in the tree T, 
its label is the absolute value of the difference between its father's and grandfather's labels (see 
Figure [SI). Notice that all labels in T are nonnegative. 

The random paths in the trees are coded by a sequence (X„)„>3 of i.i.d. random variables 
taking values in the alphabet {R, L} with probability (p, 1—p). The path starts from the root and 
goes through its only child. Then the following steps are given by {Xn)n>3- Each R corresponds 
to going through the right child (right step) and each L corresponds to going through the left 
child (left step). Note that F„ (respectively F„) is the label read in T (respectively T) at the end 
of the path X3 ■ ■ ■ X„ . 

Observe that in the linear case, each step of the path can be interpreted as the right product 
of (F„, Fn+i) by one of the two matrices 

,3, ;) o, 
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Figure 3. First lines of the trees T (left) and T (right). 

The non-linear case involves multiplication by matrices A, B and C := but their distri- 

butions is no longer i.i.d., which makes this interpretation less convenient. 

2.2. Reduction of paths. Our method relies on some properties of the trees T and T illustrated 
in Figure |4l In the following, we say an edge connecting a vertex v and its child v' is labelled by 
(a, h) when a is the label of v and h is the label of v' . 

Suppose a random path goes through an edge labelled (a, h) and then follows the pattern RLL 
in T. Then it ends with an edge which is also labelled (a, 6). Therefore, we can remove from 
{Xn)n all occurences of the pattern RLL when studying the sequence {Fn)m as long as we keep 
in mind the number of such deletions. 

The linear case is a bit more complicated. Suppose now a random path in T goes through an 
edge labelled (a, b). Notice the labels of the left and the right child are respectively b—a and b + a. 
If the path follows the pattern RLL, then it ends on a vertex labelled by —5 and whose left child's 
and right child's labels are respectively —{b + a) and —(b — a). Since we are only interested in the 
behaviour of \Fn\, this allows us to remove in (X„)„ each pattern RLL, provided we exchange the 
following letter and keep in mind the number of deletions. 

This reduction process observed in T can be translated in the language of matrices by the 
following relations satisfied by A and B: 

ABBB = -A 
ABBA = -B . 

2.3. Reduced random Fibonacci sequence in the Unear case. To formalize the reduction 
process, we associate to X3 . . . X„ a (generally) shorter word Wn = ■ ■ ■ y^^^j by the following 
induction. 

• fc(3) ^ 3 and Y;^ = X3. 

• Wn+i is deduced from Wn in two steps. 

Step 1: Add one letter (i? or L, see below) to the end of Wn- 
Step 2: If the new word ends with the suffix RLL, remove this suffix. 
Thus, we have either Wn+i = W^nyA(n)+i (^^'^ + 1) = H"^) + l)i or Wn+i = 
■ ■ ■ Ykin)^2 (and Hn + 1) = k{n) - 2). 
The letter which is added in step 1 depends on what happened when constructing Wn'- 

• If Wn was simply obtained by appending one letter to Wn-i (or if n = 1), we add Xn+i 
to the end of Wn- 
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Figure 4. The reduction pattern in the tree T (left) and in the tree T (right). 



• Otherwise, we had removed the suffix RLL when constructing Wn', we then add Xn+i to 
the end of Wn, where R := L and L := R. 

Example: Let . . .Xg = RRLLRLR. Then, the successive reduced sequence W3, . . . , Wg are 
given by R, RR, RRL, R, RL, L. 

Observe that the label read in the tree T at the end of the path coded by Wn has the same 
absolute value as Fn- 

Lemma 2.1. We denote by \Wn\R the number of R's in Wn- We have 
(4) \Wn\R >+oo a.s. 

n — ^00 

In particular, the length k{n) ofWn satisfies 

k{n) > +00. a.s. 

n — >oo 

We postpone the proof of Lemma [2. II to the end of the section. We will need in the sequel the 
following definition: 

Definition 2.2. We say that survives if, for all m > n, k{m) > k{n). 

The divergence of k{n) in the previous lemma shows that, almost surely, for any fc > 3, there 
exists a smallest Uk such that k{nk) = k and Y^*" survives. In this case, Y^'' ~ Y^ for all m > Uk. 
We will then set 

^fc ^fe • 

Note that Y3Y4 . . . contains infinitely many _R's (Lemma 12. ip . and no pattern RLL. Therefore the 
only place where consecutive i's can appear is at the beginning of the sequence. However, these 
starting L's are not relevant. If Y3Y4 . . . starts with an L, we can delete the first three letters 
without changing the values of the labels read along the path (see Figure [3]). Without loss of 
generality, we can henceforth assume that Y3 ~ R. 

Proof of Lemma \2.1\ We consider the successive changes in the number of _R's in Wn, which we 
denote by 5*1, 52, . . . G {±1}: If the j-th change corresponds to appending some R to Wn, Sj = 1; 
and if the j-th change is the deletion of some suffix RLL, Sj = ~l . Observe that if Sj = 1, the 
only way to get = —1 is to draw two L's for the following two X„'s (to remove the suffix 

RLL before another R is appended to Wn). Therefore, 

P(5,+i = = 1, . . . , 5i) = (1 - p)'. 

On the other hand, Sj = —1 means that we have just deleted some suffix RLL, so that if the next 
Xn is L (which happens with probability (1 — p)), we will have Sj+i = 1. We thus get 

TP{S,+i = l\S, = ~l,S,.i,...,Si) > (l-p). 
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We claim that the above requirements on the conditional distribution of Sj^i imply 

1 ^ 11 
(5) liminf - ^ ls,=i > > - a.s. 

This comes from the comparison between the stochastic process {Sj) and the Markov chain (A'/^) 
taking values in {±1}, satisfying M\ = 1 and having the following transition probabilities: 

P(A/,+i = -1|M, = 1) = {l-vf. 
P(A/,+i = 1|A/, = -1) = (1-p). 

Note that the invariant probability measure for the Markov chain (Af^) assigns mass 2^ to the 
point 1, so that 

1 ^ 1 
lim — \mi=\ = a.s. 

We are now going to extract a subsequence {Mg.) from [Mi) by deleting only some — I's, and 
such that {Mij) has the same distribution as (5*^). We set i\ := 1. Suppose that £j is known. If 
A/fj. ~ 1 or Mfj+i = 1, we set ^j+i Ij + 1. Otherwise, Mi. = Af^j+i = —1. Let us denote by 
(3 the probability P(S'j+i = l\Sj = Mi^,Sj-i = Af^^._,, . ..,Si= Af^J. We know that /? > 1 - p. 
We then set 

{ij + 1 with probability 
13-1 + v 
infifc > e, : Mk = 1} with probability —. 
P 

In this way, we get 

P(A/f^.^, = 1|A/,^. = -1, Me^_„. . . , Me,) = /?. 

Since the proportion of I's is greater in (Af^.) than in {Me), ([5|) follows, which in turn implies (jl]). 

□ 

2.4. Reduced random Fibonacci sequence in the non-linear case. Wc associate to . . . X„ 
the word Wn, which is obtained by the same induction as Wn, except that the letter added in 
Step 1 is always Xn+i- The label read in the tree T at the end of the path coded by Wn is equal 

to Fn- 

Lemma 2.3. Wc denote by \Wn\R the number of R's in Wn- Ifp> 1/3, wc have 

\Wn\R > +00 a.s. 

n — ^00 

In particular, the length of Wn goes to infinity almost surely 

Proof. Since each deletion of an R goes with the deletion of two L's, if p > 1/3, the law of large 
numbers ensures that the number of remaining i?'s goes to infinity. □ 

The non-linear case for p < 1/3 will be treated later (see Section [521) ■ 

2.5. Survival probability of an R. We are now able to study both cases by introducing the 
probability c of appending an R after a deletion of the pattern RLL: c ~ 1 in the linear case and 
c = p in the non-linear case. In the sequel, we consider both the linear case for any p s]0, 1] and 
the non-linear case for p > 1/3. For simplicity, wc will use the same notations Wn — ^3" . . . Yk{n} 
for the reduced word and (Yk) for the sequence of surviving letters in both cases. 

Observe that, by construction of the sequence Wn, if Y^^^ has been appended at time n, its 
survival only depends on the value of Y^^_^^ itself and the future Xn+i, Xn+2 ■ ■ ■■ We define 

Pfl := P (y^^^i survives {Y^^^^^ ~ R has been appended at time . 

A consequence of Lemma [2.11 and Lemma [2.31 is that p^ > 0. We now want to compute pu as a 
function of p. 
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We first need to analyze the following situation: Assume that in the construction of some Wng 
we have deleted a suffix RLL. Then the survival of YJ]° . depends on the nature of what we call 
the next touching letter, defined by the following algorithm: 

Step 1 Set n := ng. 

Step 2 If the letter appended to Wn is L. it may interact with Yj^^^^ to delete it. Return L as 
the next touching letter and halt the algorithm. 

Step 3 If the letter appended to Wn is R and if this R survives, then so does Y^^^ and we 
return R as the next touching letter and halt the algorithm. 

Step 4 Else, the letter appended to Wn is a non-surviving R. Then there exists a smallest 
integer m > n such that k{m) = k{n) (corresponding to the time when this R is deleted). Then 
set n := m and go back to Step 2. 

Each time the algorithm enters Step 2, it has a probability 1 — c to directly return L. Since 
1 — c > 0, the algorithm ultimately halts with probability 1. 

We now go back to our computation of p^. Here is an exhaustive list of all cases in which 
^k{n) ~ a-ppcnded at time n, survives: 

• Case 1: R [41^] R 

'^ki'n+i) ~ "'^n+i = Either it survives or it does not survive but the first touching letter 
after £ > 1 deletions is an R\ This happens with probability 

(6) Pi Wi? +P ^(1 -Pi?)^c^Pi? 



l-c{l-pn) 



Case 2: R [R^] L [R^] R 

Either Y^^^^^ = Xn+i ~ L, or Y^,^]^^^ = Xn+i = R which does not survive but the first 
touching letter after (.>! deletions is an L. In both cases, this L is immediately followed 
either by a surviving i?, or by a non-surviving R and the second touching letter is an R] 
In view of the probability computed in Case 1 , this happens with probability 

P2:= \{l-p)+pY.{l-pnYc'-\l-c)^ 



Writing pji = pi + P2, we get that pn satisfies 



1-c{1-pr) 



Pr{c^Pr + + 2c(l -p- c))pR + (1 - c - 2p)(l - c)) = 0. 

When pfi 7^ 0, which is true in the linear case for any p and in the non-linear case for p > 1/3, 
this equation has only one non-negative solution given by 



-p2 - 2c(l -p-c)+ py/p^ + 4c(l - p) 

(7) PR^ . 

Observe that in the non-linear case for p < 1/3, the above expression is non-positive. Thus. 
PR = 0. 

2.6. Distribution of the reduced sequence. We deduce from the preceding analysis the prob- 
ability distribution, when p^ > 0, of the reduced sequence Y which is the concatenation of all 
surviving letters. By construction, since we assumed Y^ ~ R, consecutive Us are not allowed in 
Y, hence 

p(yfc+i = R\Yk = i) = 1. 

Moreover, if = R has been appended at time n and survives, the following letter yA;(ri)+i in 
Y only depends on Xn+i, X„+2, ■ • ■, which implies that (Yk) is a Markov chain. Observe that if 
Yfe = R, it is followed in Y by another R only in the last two cases of our exhaustive list. Hence, 

(8) a:=I'{Yk+i^R\Y,=R) 



PR 2c-p+ y/p'^ -I- 4c(l - p) 
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which leads to 



(9) 



3p - 2 + V5p2 -Sp + A 

(hnear case , 

2p ^ ' 

^ (non-hnear case). 



.P+\/p(4-3p) 

The invariant probabihty measure of this Markov chain is given by 

1 I- a 



(10) (Aifl,ML) 



2-a' 2-a 



^From now on, we denote by the probabihty distribution on {R,L}^^'^'^^^^ under which (Yfe) 
is a Markov chain with = i?, 

P„(yfc+i = R\Yk - i?) = a and P„(yfc+i = R\Yk = L) = I . 

2.7. Compression rate. We are also interested in the ratio of surviving letters when k is a 
large integer. The number Sk of i?'s in Y3 . . . Yk, which is the number of surviving i?'s up to time 
Uk, satisfies Sk/k = fin + o(l). Let dk be the number of deleted R's up to time Uk- Observe that 
the total number of R's drawn up to time Uk is Sk + dk, so that Sk/ {sk + dk) ~ Pr + o(l). Since 
each deletion of an R comes with the deletion of two L's, we have Uk = k + 3dk- Therefore, 

(11) .a:^(l + 3^in^^^y' a.s. 

Uk fc^oo \ Pj^ J 

P 

It will be useful to see cr as a function of a. ^From ([6|), we get a = pi/pR = -. r. Hence, 



l-c(l-pfl)- 



1 p 

Pr = ■ 



a 



Moreover, ([5]) yields 

il-af 



p = a ~ c- 



a 



Taking (|10p into account, an elementary computation leads to 

(2a-l)(2-a) 

(1^) '^^ — 2 n — 

3. Continued fractions in the tree R 

3.1. The tree R. The reduction of the sequence (X„) lead us, both in the linear and the non- 
linear case, to the study of a Markov chain (Yk) whose distribution is P^ on {R, L}^^'^ - -^ . The 
only difference between the linear and the non-linear case is the value of the parameter a. 

Let us consider the random sequence of integers (Gfe)fe>i, where Gi = G2 = 1 and Gk is the 
label read in the tree T (T in the non-linear case) when following the path coded by Y3 • • • . 
We get that for all fc > 1, Gk has the same absolute value as Fn^. {Fn^. in the non- linear case). 
We are thus left with the estimation of the exponential growth of the reduced Fibonacci sequence 
{Gk)k>i- 

Since (Yk)k has no pattern RLL, the reduced Fibonacci sequence (Gfc)fe>i can be read along a 
random path in the tree R introduced in [8], which is a sub-tree of T and T. The tree R is defined 
as follows: The root of R has only one right child, which itself has only one right child. Any other 
vertex v has either one right child or two (left and right) children, depending on whether v is itself 
a left child or not: A left child has only one child, whereas a right child has two children. Each 
vertex is labelled with an integer : The root and its only child are labelled by 1. If w is a right 
child, its label is the sum of its father's and grandfather's labels; If w is a left child, its label is the 
difference between its father's and grandfather's labels. (See Figure O) 

Please note that the step from the root to its child does not appear in Y: The condition Y3 = i? 
corresponds to the fact that the only child of the root has only one right child. For technical 
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2 4 4 2 8 

\ V\ /\ \ /\ 

5 1 7 3 5 7 3 13 

Figure 5. First lines of the tree R. The random path marked with the bold 
edges is coded by the sequence Y = RLRRL . . . 




Figure 6. Distribution probability of the random path in R. 



reasons, we will sometimes need to add an extra R at the beginning of Y , representing the step 
from the root to its child, but this will be done explicitely. 

From the preceding section, we know that the distribution of the random path is given by a 
Markov chain: Each left step is followed by a right step, and each right step is followed by a right 
step with probability a and by a left step with probability 1 — a, where a is given by ([8|) (see 
Figure O. 

Let us recall some important properties of the tree R. First, it is easily proved by induction 
that the label of each vertex is a positive integer. As a consequence, we get that 

Vfc > 1, Gfe > 0. 

Another easy induction shows that, if x and y are the labels of a child and its father, then x and 
y are relatively prime. Moreover, for each pair (a;,?/) of relatively prime integers, there is exactly 
one vertex in R which is labelled by x and its father by y (see [5]). Therefore, for any positive 
rational number g, there is a unique vertex a; in R such that q is the quotient of the label of x by 
its father's. This vertex a; is a right child if and only li q> 1. 



For all k > 3, let Qk ■— Gk/Gk-i, so that 



(13) 



^logGfc = T^^ogQi 



i=3 
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The exponential growth of Gk will thus be deduced from the probability distribution of {Qk)k>3-, 
which is related to the development in continued fractions of Qk- 

3.2. Shape of a path and continued fractions. Let ^ be the set of all finite sequences 
y = VaVi ■ ■ - Uk (k > 3) of R's and L's with = R and no pattern LL. Each sequence in 3^ can 
be interpreted (with the same conventions as above) as a finite-length path in the tree R. 

To each sequence y — y^y^ . . .yk & we associate a rational number q{y) defined as follows. 
We decompose the path ysy^ ■ ■ .yk into pieces which are either elbows RL or single right steps R. 
(This can be done in a unique way for all y G ^.) Next, we introduce a cutting between each 
pair of (successive) identical pieces. We thus obtain a partition of the path into £ blocks; Let ai 
be the number of pieces in the last block, 02 the number of pieces in the last but one block, and 
so on. If the last piece of the last block is an elbow, q(y) is given by its development in continued 
fractions 

qiy) := [0, ai, . . . , Of] = + — < 1, 



1 



ai 



1 



02 



otherwise, q{y) is set to 

q{y) := [oi, . . . , a^] = ai H > 1. 



a2 



1 



The following proposition shows that the random variables (Qk) are precisely given by the 
preceding computation. 

Proposition 3.1. Let y = j/3 . . . j/fe G coding a Gnite-lcngth path in R. Let gi — gi — 

1, ■ ■ ■ , gk be the labels of the vertices visited by the path. Then q{Ry) = gk/gk-i- 

(See Figure El) 

Proof. We proceed by induction on k. Let us denote by qk the quotient gk/gk-i- For fc = 3, the 
sequence Ry is reduced to RR: We thus have two single right steps separated by a cutting and 
we check that 

q3='^ = l + j=qiRR). 

Assume the result is proved up to fc — 1 and consider a sequence y oi R and L's of length k. Let 
us analyze the three possible configurations for yk-iyk- 

• yk-iUk = LR (see Figure HI case (1)): Since the last letter is an R, we have 

5fe gk-i + gk-2 , , 1 
qk = = = 1 



5fe-l 9k-l qk-l 

On the other hand, the induction hypothesis gives qk-i ~ q^Rya ■ ■ ■ Uk-i) ~ [0, ai, 02, . . . , ae\ 
with [ai) given by the number of pieces in the blocks. Hence, q^. = [ai + l,a2, . . . ,a^]. 
Since there is no change in the cuttings when appending the last R, the number ai of 
pieces in the last block is increased by 1 and the other ones are left unchanged. Therefore, 
q{Ryz ■■■yk) = qk- 

Uk-iVk = RR (see Figure[8l case (2)): We still have gfc = 1 + 1/qk-i- The difference with 
the previous case is that qk-i = [0.1,0,2, ... ,ae] > 1, thus qk ~ [1, ai, 02, . . . , a^]. Here, the 
last R has introduced a new cutting, so that we have one more block of length 1 and the 
other blocks are left unchanged. Therefore, q{Rya . . .yk) = qk- 
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04 = 6 



ag = 1 



as = 2 



ai = 1 



Figure 7. Consider the sequence y = RLRRLRRLRLRLRR £ ^. The corre- 
sponding path in the tree R is represented above, as wch as {gj)i<j<i6- After 
decomposing Ry into elbows and single right steps, and introducing the cuttings 
between identical pieces, we obtain the development in continued fractions of 
q{Ry) = 515/515 — 27/20 by counting the number of pieces in each block: 27/20 = 
[1,2,1,6]. 



yk-Wk = RL (see Figure |51 cases (3) and (4)): Since the last letter is an L, we have 



Qk 



9k 
Qk-i 



9k-l — 9k-2 
9k~l 



1 



qk-i - 1 



The induction hypothesis gives that qk~i = qiRys ■ ■ ■ Uk-i) = [o-ij 02, . . . , ag], with (a^) 
given by the number of pieces in the blocks. Hence, the development in continued fractions 
of Qk depends on the value of ai. 

If ai = 1 (case (3)), we get qk = [0, 02 + 1, as, . . . , a^]; Appending the last L transforms 
the last piece into an elbow. The fact that ai = 1 means we had a cutting just before this 
last piece, which disappears after the transformation. Therefore, we have one less block 
and the number of pieces of the last remaining block is increased by 1. 

If ai > 1 (case (4)), we get qk = [0, 1, ai — 1, 02, ... , ag]; The fact that ai > 1 means we 
had no cutting just before the last piece, and one is created when appending the last L. 
Therefore, we have one more block of length 1, and the number of pieces of the last but 
one block is decreased by 1. 

In both cases, we get q{Ry3 ■ ■ - yk) ~ qk- 

□ 
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(1) (2) (3) (4) 

Figure 8. From qk-i to qk- (1) corresponds to the case when qk-i < 1, (2), (3) 
and (4) correspond to the cases when qk^i > 1- 

3.3. Stern-Brocot intervals. We recall here the definition of the positive Stern-Brocot inter- 
vals, which are the real intervals of the form [j, ^] where a,b,c,d are nonnegative integers with 
ad — be = —1. (We adopt the usual convention that ^ = +cx).) All these intervals can be obtained 
inductively, starting with the interval [j, i] (the only Stern-Brocot interval of rank 0). Each Stern- 
Brocot interval [|, ^] of rank r is cut into two Stern-Brocot intervals of rank (r + 1): [|, and 
[f+5' l}' -f^ccall also that is called the mediant of | and ^. 

Let us fix a sequence s = si . . . Sm G ^ (we also include here the case where s is the empty 
sequence) . We now consider the set of finite paths in the tree R having s as a suffix: 

{Ry : y = {yi)3<i<k e ^, k>m + 2; yk-m+i . . . i/fc = si . . . s^} . 

Proposition 3.2. Let r be the number of pieces in the decomposition of s into elbows RL and 
single right steps R. Then 

• The closure q{&'s) of q{'3^s) is a Stern-Brocot interval of rank r; 

• if r is even (respectively odd), q{'3^RLs) is the left (respectively right) Stern-Brocot sub- 
interval of rank r + 1 of q{'3^s), s^nd q{'3^Rs) is the right (respectively left) one. 

Proof. We proceed by induction on r. If r = 0, s is the empty sequence and is the set of all 
finite paths in R, so (?(^%) = [0, +oo]. Moreover, ^Vrs = % is the set of all finite paths ending 
with an R, hence q{'3^Rs) = [l,+oo]. 'S^rls = ^rl is the set of all finite paths ending with an 
elbow and q{^RLs) ~ [0, 1]. Therefore q{?^RLs) is the left Stern-Brocot sub-interval of rank 1. 

We consider a suffix s such that qi^K) is a Stern-Brocot interval of rank r, and its two possible 
extensions Rs and RLs. Assume for simplicity that s ends with an R (the proof is the same when 
s ends with RL). Let £ be the number of blocks in s: q{s) — [ai, . . . , a^]. Since the last letters of 
y give the beginning of the development in continued fractions of q{y), all q(y) for y € have 
their first {£ — 1) partial quotients fixed, and equal to those of q{s). Moreover, their £-th partial 
quotient is at least at. Thus, q{'3^s) is the interval whose bounds are [ai, . . . , a£_i] and [ai, . . . , a^]. 
Notice the way these bounds are ordered depends on the parity of t. [ai, . . . , a^-i] < [ai, . . . ,ai] 
iff £ is even. 

Either Rs (if the first piece of s is an R) or RLs (if the first piece of s is an RL) counts 
one more cutting than s and gives rise to the interval whose bounds are [oi, . . . , a^-i, ai, 1] and 
[ai, . . . ,ai-i,ai\. The other one gives rise to the interval whose bounds are [ai, . . . ,af_i] and 
[oi, . . . , ag-i, a£ + l]. In fact, since [ai, . . . , a£_i, a^, 1] = [ai, . . . , a£_i, a^-f-l], the two intervals have 
a common bound. Wc let the reader check that this common bound is the mediant of [ai, . . . , a^_i] 
and [ai, . . . ,ae], so that q{'3^Rs) and q{'3^RLs) are two Stern-Brocot intervals of rank r As a 
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consequence of the previous remarks, the table below gives the relative positions of q{'3^RLs) and 





the first piece of s is i? 


the first piece of s is RL 


i odd 


Rs ^ left 
RLs ^ right 


RLs ^ left 
Rs ^ right 


I even 


RLs ^ left 
Rs ^ right 


Rs ^ left 
RLs ^ right 



Each time we extend the suffix s with one more piece, either this piece is similar to the first piece 
of s, which introduces a new cutting, and £ is increased by 1, or the new piece is different and £ is 
unchanged. Therefore, we move either vertically or horizontally in the table, so that the relative 
positions of q{'3^Rs) and q{'3^RLs) alternate. □ 



3.4. Probability distribution of {Qk)k>3- We now turn back to the Markov chain Y following 
the probability distribution P„ (see Section [276| . The ergodic theorem for this Markov chain gives 
that, almost surely, 

k 



RY3...Y,e?yn ^ > Aifl- 

3 



If we decompose the sequence Y3 . . .Yi into pieces RL and R, it is not hard to see that all but the 
last piece appear independently, with probability 1 — a for RL and a for R. Therefore, if we fix 
some s & ^ and denote by \s\rl (respectively |s|_r) the number of pieces RL (respectively R) in 
the decomposition of s into pieces, we get from the law of large numbers that 



1 

(14) -V W3...r.G;% >c(s)al-'l«(l-a)l 

K ^ — ' k^oo 
i=3 

where 

cis) 



l-tR/a if s ends with an R, 

(1 — fiR)/{l — a) otherwise. 



Since for all i > 3 we have Qi ~ q{RY^ ... y^), it is natural to introduce the following probability 
distribution ly^ on ]R+: i>a is defined by 



In view of Proposition 13. 2i this amounts to define it inductively on Stern-Brocot intervals in the 
following way: First assign mass 1 — a to [0,1] and a to [l,oo]. Once Va is defined on some 
Stern-Brocot interval of rank r, a proportion a of its mass is given to the left (respectively right) 
subinterval of rank r + 1 when r is odd (respectively even) (See Figured]). We can notice the 
similarity between this construction and the Denjoy-Minkowski measure /i^"-* presented in [2]. 
The difference lies in the fact that the proportion a is in our case alternatively given to the left 
and the right subinterval. 

From ([T4|). we obtain that, for all / e L^(z^q), 



-1 ^ -1 fl fOC 

(15) 7:E/(Q^)^ .f{x)du^{x) + ^ f{x)du^[x) 

k ^ fc^oo 1- a Jq a Ji 



a.s. 



Remark 3.3. Observe that we need the correction c{s) in p4)) because the last piece of Y3 ... 1^; 
has a different distribution. If we first cut the infinite sequence (li)i>3 into pieces R and RL, 
and denote by ij the index at the end of the j-th piece, then the distribution of q{RY3 . . .Yi^) 
converges to Va without correction. 
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4. The Lyapunov exponent 



4.1. Computation of the Lyapunov exponent. We now end the proof of Theorem ll.il It is 

easily seen tliat x i-^ log a; belongs to L^{i'a)- Using ^3]) and ([15]), we obtain 



-logGfc » 



Using (fTTj) and fin 



\ — a 
1/(2 — a), we get 



log(a;) dva{x) + 



(16) 



1 



log (a;) dvaix) 



a.s. 



log(x) (il/a(x) H / log{x) dl^aix) 



In the linear case, since we are dealing with a product of i.i.d. matrices, we know that the limit 
of {1/n) log \ Fn\ exists almost surely, and is given by the largest Lyapunov exponent. 

Of course, we get the same formula as ((T6| for (l/iik) logi^„j. in the non-linear case for p > 1/3 
(where a is given by ([9])). As we already pointed out, this case cannot be reduced to a product of 
i.i.d. matrices, therefore we need a little argument to get the almost-sure existence of the limit of 
(1/71) logF„. It consists in controlling the size of the deleted blocks. 

We already know the almost-sure convergence of (l/rifc) logi^„j^, along the subsequence (uk) 
corresponding to surviving letters after the reduction process. Consider now n lying between Uk 
and rifc+i. We have 



-logFn - — log Fn, 
n Uk 



< 



Uk 1 



nk 



log 



F„. 



We need to control the quantity Tk :~ — nk, which is 1 plus the number of deleted letters 
between two successive surviving letters. The probability distribution of the random variable Tk 
is given by the law of the following stopping time for the i.i.d. sequence on the alphabet {i?, L} 
with probability {p, I — p): Draw a sample of this process and stop the first time the number of 
L's is equal to twice the number of R's plus one. This sample without the last L corresponds 
to all the patterns RLL, between two successive surviving letters, which are removed during the 
reduction process. When p > 1/3, the stopping time is almost surely finite and its expected value 
is finite. 

Since the T/c's are i.i.d. and have a finite expected value, we get 

n-Hk ^ Tk 

k k- 



< > a.s. 



n 



Observing that F„ < maxji^n^ , F„j,_j }2" the convergence along the subsequence (n^) is 
enough to conclude that the limit of (l/ri)logF„ exists almost surely, and is given by the right- 
handside of (1161) . 

To this end, we use some 



It remains now to prove that this limit is equal to log (a;) dva 
changes of variables in the computation of the integral. 



Lemma 4.1. For all f E L^{va), 



1 



/ 



1 -a; 
1 -a; 

X 



dVa{x) 
dVa{x) 



{I -a) 



f{x)dVa{x) 

X) 

f{x)dVa{x) 



Proof. It is sufficient to prove the equalities when / = l[a/b.c/d] is the indicator function of a 
Stern-Brocot interval {ad — he = —1). The first equality becomes, for a>h, 



1 - a 



a c 
b'd 
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Observe that [^^, ^-^] is also a Stern-Brocot interval, which has the same rank r as [f, §]• For 
r = 1, the equality holds trivially. Since x ^ 1/(1 — x) is increasing, the right sub-interval of [|, ^] 
corresponds to the right sub-interval of [^^, ^7^] . The result follows by induction on r. 
For / = Ir. 



the second equality reduces to 
d b 



c + d^ a + b 



If ^] is a Stern-Brocot interval of rank r. 



(1 - a) Va 



a c 
b' d 



is also a Stern-Brocot interval, but of rank 
1. For r = 1, the equality holds trivially. Since x {1 — x)/x is decreasing, the right sub- 



c+d ' a+b 



interval of [|-, ^] corresponds to the left sub-interval of 
on r. 



d b 
c+d' a+b 



The result follows by induction 

□ 



Applying Lemma 14.11 to x log x and summing, we get that 

\og(x)diyaix) ^ / log{x)diya{x) + (l — a) / log{x)di^aix)- 

01 Jl Jo 

It follows immediately that 

log(a::)dz^Q(a;) = / \og{x) di^ a {x), 



(17) 



1 - 2a 7o 

log(a;)dz^Q(a;) = / log(x)di^c<(a^)- 

1 1 - 2q; Jo 



Substituting in the right-handside of and recalling p^ . we conclude the proof of Theorem ll.il 
for all p in the linear case and for p > 1/3 in the non-linear case. The non-linear case for p < 1/3 
is treated in the next section. 

4.2. Variation properties of the Lyapunov exponent. We now want to prove that p 1— *■ 7p 
(non-linear case) is a non-decreasing function. We first establish a comparison lemma. 

Lemma 4.2. Let x be a path in the tree T, and let x' be obtained from x by turning an L into 
an R. Then any label read along x is always smaller than the corresponding label read along x' . 

Proof. Let y ~ . . .y^ G coding the end of a finite-length path in R. Assume that the 
vertices of the edge preceding y are labelled by a and b. Then the last vertex of the path is 
labelled by a linear combination of a and b of the form d{y)a -\- n{y)b, where n(i/) and d{y) are 
integers depending only on y (see Figure IHl left). More formally, n{y) and d{y) can be defined by 
the following induction. n(0) := 1, d(0) :== 0, n{R) d{R) := 1, and for i > 3 



n(y3 ■ 
d{yz ■ 
n{y3 ■ 

d{y3 ■ 



1, d(0) 

■ v^R) 

■ y.L) 

■ y.L) 



n{y3 ■ 
d{y3 ■ 
n{y3 ■ 
d{y3 ■ 



■ yi) + n{y3 . 
■yi) + d{y3 . . 
■yi) - n{y3 ■ 
■yi) - d{y3 . . 



■ yt-i) 
■yi~i) 
■yt-i) 

■ yi-i) ■ 



Since y codes a path in R, an induction shows that d{y) and n{y) are nonnegative. 

Consider two paths x and x' in T differing at level i : We decompose the end of the paths from 
level i as LU'y and RL^y, where y starts with an R and r > 0. 

Suppose first that, after the difference, all letters of x and x' are i's (y = 0). We let the reader 
check that the labels after level i and i + 1 are well-ordered. Moreover, the label after level k > i + 2 
in x' is equal to the label after level fc — 3 in x. This can be seen by making a reduction (removing 
a pattern RLL) in the path a;'. Denoting by a and b the (nonnegative) labels after level k — 3 and 
A: — 2 in a;, the label after level fc in a: is given by |5 — |5 — a|| < a. 

Suppose now that the sufhx y is reduced (y G ^). The above argument shows that labels 
after level i + r — 1 and i + r are well-ordered: Denote these labels by a and 6 in x, and a' and 
x' in x' . Then the label at the end of x (respectively x') is the linear combination d{y)a + n{y)b 
(respectively d{y)a' + n{y)b'). Since d{y) and n{y) are nonnegative, we conclude this case. 
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h d(0) = 0, n(0) = 1 

^+6 d{R) = l,n{R) = l 



d{y)a + n{y)b 

Figure 9. Definition of d{y) and n{y) (left); Comparison between two paths in R (right). 



In the general case, we make all possible reductions on y. We are left either with a sequence in 
^ or with a sequence of Us, which are the two situations we have already studied. 

□ 

Proposition 4.3. Ifp< 1/3, jp :— lim„(l/n) logF„ exists and is equal to 0. The function p^"fp 
is continuous and non- decreasing. 

Proof. We start by proving that Vnnp\^i/^^p = 0. We can extend the definition of i^a to any 
a e [0, 1], and also set for < a < 1/2 



7(a) 



Observe that, if [a,b] is a Stern-Brocot interval, then [1/fe, 1/a] is also a Stern- Brocot interval, 
satisfying 

i^a([a,b]j = i^i_a(^[l/6, 1/a]^. 

Therefore, we get 

7(a) = -7(1 - a). 

In particular, 7(1/2) = 0. Moreover, a 1— > 7(a) is easily seen to be a continuous function (we can 
write it as a uniform limit of continuous functions). Hence 

lim 7„ = lim 7(a) = 0. 

p\l/3 

Now, let < j5 < p' < 1. Let (X„) and (^4) be random paths in T for the respective parameters 
p andp'. We can realize a coupling of (X„) and {X^) such that for any n, Xn = R implies X!^ = R. 
From Lemma [4.21 it follows that the label Fn read along X is always smaller than the label 
read along X'. If we choose p < 1/3 and p' > 1/3, we get that 



lim sup — log F„ < lim — log 

n n 



Since limp/\^i/3 7p' = 0, we deduce that 7^ = for any p < 1/3. Moreover, this argument obviously 
shows that p i— > 7^ is a non-decreasing function. □ 

Corollary 4.4. The function p^ ^p is increasing on ]0, 1[. Tiie function p^^p is increasing on 
]1/3,1]. 



Proof. Recall that 7(a) = \ogxdva{x) = 7p, for a = 2p/{p+ ■\/p(4 — 3p)). The function 
p i-^ 2p / {p+ p{A: — ip)) is increasing and sends ]l/3, 1] onto ]l/2, 1]. Hence 7(a) is non-decreasing 
on]l/2,l]. 

In the linear case, we also have 7p = 7(a) where now a = (3p— 2-|--\/5p^ — Sp + 4)/(2p) S]l/2, 1]. 
Since this expression is also increasing va.p, p^ jp is non-decreasing. Moreover, we easily deduce 
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from pT)) that jp > for any p e]0, 1[. Wc also know from 7. that 7^ is an analytic function of 
p G]0, 1[, thus it is increasing. 

This in turn implies that 7(0;) is increasing on ]l/2, 1], and so is jp for p g]1/3, 1]. □ 

4.3. Derivative of the Lyapunov exponents. The following proposition gives a formula for 
the derivative of 7 with respect to a, which uses the product measure Va'Si J^a- 



Proposition 4.5. For all 1/2 < a < 1, 

(18) 7 (a) = 7(a)- 



2a - 1 



1 



log 



X + y + xy 
x + y + 1 



Va ® Va{dx, dy) 



(2a-l)(a2-a + l) - 

Before proving this formula, we now use it to compute the derivatives of the Lyapunov exponents 
for p = 1. When p = 1, a = 1 (both in the linear and in the non- linear cases), and is the Dirac 
measure on (j). (fTS]) yields 

logs 



-(l) = log0 + log^^ 



da , 



Since (1) = 1 (both in the linear and in the non-linear cases), we get 
dp 

^(l) = ^(l) = (log5)/2. 
dp dp 

We now turn to the proof of Proposition [131 

Proof of Proposition \4 -SI We fixl/2<Q; — e<a<l. Let F be a Markov chain following Fq, 
decomposed into pieces R and RL. We decide independently to change each piece R into RL 
with probability e/a. We thus obtain a new Markov chain Y' following IPa-e- Let {Gk)k>i and 
(G"j,)fc>i be the labels read along the paths Y and Y' respectively. We introduce the subsequences 
(km) and {k'^) such that Gk^ and GJ,, are the labels read after the m-th piece of Y and Y' 
respectively. We know that 



7(a)- 7(0; -e) 

(19) 
Since 



(j{a) lim E 



T-l0gGfc„ 



— a{a — e) lim E 



{(j{a) ~ a{a — £)) ^\ \ + ^(a ■ 



cr(a) 



e) lim E 



T^logGfc„ 



TT^ogG'k, 



lim E 



h 

^ >2 

m m~fOO 

the last term in (fT9)) becomes 

a{a — e) 

2 — a m-*co 

Therefore, we obtain 

7(a) — 7(0; - 



(20) 

An easy computation shows that 

' o{a) — a{a — e) 7(0;) 



a and ^ 



m m^oo 



— logT^ 



2 — a + e a.s., 



7(a-e). 



2- a 



u{a) — a{a — e) "f{a) 

a(a~-e) 1 ^ 
+ 7^ T - lim E 

(Z — a) £ m^oo 



-7(a - e) 



(21) 



lim 



<7{a) 2 — a 



7(a-e) =7(a) 



1 + 2a - 2a2 
'(2a- l)(a2 - a + 1) 



Let us now turn to the last term on the right handside of (PD)) . For all i > 1, let us consider 
the «-th piece which is different in Y and Y' . We denote by rrii the number of pieces which have 
been seen before the i-th change. Let at := Gk^.-i and bi := Gk^. be the labels along Y of the 
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vertices of the edge preceding the i-th change. Similarly, let := Gk' 
corresponding labels in Y'. Using the fact that i/rrii —* e a.s., we get 



.1 and b', G',, 



be the 



(22) 



lim E 

m — >-oo 



1 1 Gk„ 



: hm E 

■i— *CXj 



1 , fli 
-log — 



Between the i-th change and the edge (0^+1,6^+1), Y and Y' share a common part (see 
Figure fTO|) . We set rii := n(y*^'+^) and di := (i(y^'+^), where n(-) and d(-) are the notations 
introduced in the proof of Lemma We then have the following induction. 



This yields 
E 



1 , ai 
-log — 



djbj + ni{a.i + hj) 



bj n.j ( bi 
- + ^ 1 + - 



1 



E 



log(^ 

a,- 



E 



log 1 



6,, 



fli + bi 



c^ + K 



Y 



a,+i = dibi + ni{ai + b^) 



&i+l 



a-+i =dj(a^ + 5-) + '^*a» 



Figure 10. The paths Y and Y' between the i-th and {i + l)-th change. 



Observe that bj/aj and rij/dj are independent. For all j, Uj/dj has a probability distribution 
r/^ which only depends on a and e. Moreover, we know that bj/aj converges in law to Va (see 
Remark 13.31). It follows that 



E 



log 



Similarly, 



E 



ttj dj 



log 1 4- 



1 



We thus obtain 



E 



1 , ai 
-log — 



logix + y {1 + x)) diyaix)dv^{y) . 



log (1 + a- + y) dva{x)dvl{y) 



X + y{l+x) / N , ^ 

log dva{x)dvl{y) . 

\ + x-\- y 
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The probability ly^ is the distribution of n{Y^~*'^)/d{Y^^'^). When £ ^ 0, the length of the 
common part yi~*2 gQ^g infinity almost surely. 

Lemma 4.6. Let Y follow the probability distribution Pq. Then n{Yi . . .Yk)/d{Yi . . .Yk) has 
almost surely a limit as k ^ oo, which follows the probability distribution Va- 

This lemma ensures that the preceding integral goes, as e — > 0, to 

x + y{l + x) 

log—— ■ dva{x)dva{y) ■ 

I + X + y 

Together with ([21]) and this achieves the proof. □ 



Y, 



da + nb 
d'a + n'b 



^{d + d')a + {n + n')b 
da + nb 

Figure 11. Labels on the path Y. 

Proof of Lemma \4^ We decompose the path Y into pieces R or RL. Recall that, under Pq,, each 
piece appears independently with probability a for 7? and 1 — a for RL. To each piece, we associate 
a real interval: Suppose the piece ends in Yk . Then the bounds of the corresponding interval are de- 
fined as n(Y3 . . . Yk-i)/d{Y3 . . . Ife-i) and n{Y3 . . . Yk)/d{Y3 . . . Yk). Observe that the first interval 
is [0, 1] with probability 1 — a and [1, oo] with probability a. If n/d ~ niY^ . . . Yk-i)/d{Y3 . . . Yfe_i) 
and n'/d' = n(l3 . . . Yk)/d{Y3 . . . Yk) are the bounds of the interval associated to the j-th piece, 
then the bounds of the interval associated to the (j-l-l)-th piece are either n'/d' and {n+n')/ {d+d') 
with probability a, or (n + n')/{d + d') and n/d with probability 1 — a (see Figure [TT]). There- 
fore, we get a decreasing sequence of Stern-Brocot intervals converging to a point following the 
probability distribution Va. O 

5. Link between i^a and Furstenberg's invariant measure 
In the linear case. Furstenberg's formula gives 

1 , , , f f , \\xA\\ , \\xB\\\ , , , 

- log \F^\ > 7p = / Ulog \^ + 1 - P log \^ dvi{x) , 

n n^oo J \ \\x\\ \\x\\ ) 

where A and B arc the matrices given in ^ and J^/ is the invariant measure on the set P(P^) of 
directions in the plane for the random walk that sends x to 'xA with probability -p and to 'xB with 
probability 1 — p. (In the above formula, x stands for any nonzero vector with direction x.) 

Directions x can be parametrized using slopes a; = (1, m) with m G (— oo, oo]. Therefore, in this 
context, V j is the probability distribution on (— oo, oo] such that, for any non-negative measurable 
function 

/,9d^/ = 1 \P9 (l + ^) + (1-^^)5 - ^) } • 

Observe that, in view of the particular form of the matrices A and B and the fact that is 
invariant, Furstenberg's formula reduces to 

Iv = / log!™! «i^/("^) + / plog— ^— f h(l-p)log— -— 2i_Ml/y(TO) 



|(l,m)l| ^ ^ 11(1, "^^ll 



+ 00 

log |m| diyf{m) . 
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(It is worth remarking that this simphfication is always vahd when deahng with hnear recurrence 
equations.) 

The difficult part is to identify the invariant measure Vf. This was done by Viswanath in the 



case p ~ 1/2, and we could note that ;//(. |IR^ 



for a = a(l/2) = (j)- 



This observation, 

together with the fact that the above equation looks very similar to our formula for 7^, made us 
suspect a relationship between Va and Vf. 

We came to the following heuristics: In the tree T, all edges whose labels have opposite signs 
correspond to a step which will appear at the end of a deleted pattern RLL. Therefore, we expect 

1-a (l-a)2 



1^/(11 ) to be equal to the frequency of deletions 



Moreover, all edges whose 



labels have same signs can be seen as belonging to a tree R. Therefore, we expect i^f{ . |R^) to 
be equal to (as in the case p = 1/2). The invariance property of i^/ for the indicator function 
of a Stcrn-Brocot interval [a, b] with 1 < a < 5 now yields 



iyfi[a,b]) =pvf 



1 a - 1 



+ {l-p)iyf 



I'fe-l 



Observe that 



1 1 

b-l ' a~l 



and 



-1 

b-l 



are also Stern-Brocot intervals, the latter lying in ]R . 

This equation is thus enough to get the measure of all Stern-Brocot intervals in ]R~. We obtain 
the measure described in Figure [T2l 



^^^-3/1 -2/1-3/2 ^1/1-2/3 -l/2_i/3°^^ 1/3 2/3 3/2 2/1 3,1 1/0 



m a(l — a) 



m+(l - a)2 



m"'"a(l — a\ 



m a(l — a) 



m {1 — a)^ 



TO (1 — a) 



m'^a{l — a) 



TO+(1 — a) 



m+ = 1 — m 



Figure 12. The measure ff on Stern-Brocot intervals of rank 0, 1, 2. First assign 



mass TO 



(l-«)2 



to 



, 0] and 



— to [0,00]. Once i>a 



— a + 1 ' ' a" ~ a. 

is defined on some Stern-Brocot interval of rank r in 1R+, a proportion a. of its 

mass is given to the left (respectively right) subinterval of rank r -f 1 when r is 
odd (respectively even). In R~, exchange a and 1 — a. 



Conversely, we easily check that this measure satisfies the invariance property. 

6. Further developments and open questions 

6.1. Extension to Viswanath's setting. The sequence (i„) studied by Viswanath in [TU] is 
defined by ti := t2 1, and 

in = ±tn-\ ± tn-2 (w > 2), 

where each ± sign is independent and takes value + or — with probability 1/2. A natural extension 
of this setting would be to choose the signs independently both with probability p for -I- and 1 — p 
for — . This definition is equivalent to ours only in the case p = 1/2, which explains why the graph 
of the Lyapunov exponent drawn on Figure 5 in |10j is different from our Figure [2l In Viswanath's 
setting, no explicit formula is known to compute the Lyapunov exponent. Can our method be 
extended to this setting ? 



HOW DO RANDOM FIBONACCI SEQUENCES GROW? 



21 



6.2. Random Fibonacci sequences with multiplicative coefficient. Consider the gener- 
alization of the random Fibonacci sequence (i^„)„ defined by -F!„+i = XFn ± F„-i for a fixed 
parameter A. 

For special values of this parameter, namely Xk — 2cos(7r/fc), fc > 3, we hope to generalize our 
method. (Observe that the present study corresponds to k = 3). We expect the reduced sequences 
to be obtained by removing patterns RL'^'^, which should give rise to the study of a Markov chain 
of order k — 2. The correspondance between random Fibonacci sequences and continued fractions 
extends to these new sequences by considering Rosen continued fractions, introduced by Rosen 
in [9] . A Rosen continued fraction expansion of a real number a; is a continued fraction in which 
partial quotients belong to AfcZ* (instead of Z!J^ as in the classical case). The values are known 
to be the only ones for which the corresponding Mobius group generated by the transformations 
z I — > z — Xk and z i — > 1 /z acts discontinuously on H O [6] . 

6.3. Average growth rate. We announced in Remark 11.21 that we were able to give an explicit 

expression for the limit of — log(E(F„)). Strangely enough, the similar question in the linear case 

n 

seems to be more difficult. We do not know yet whether the combinatorial methods used in [8] 
can be extended to this case. 

6.4. Variation properties. We could expect the formula ([T5)) to give the value of the derivative 
when p — !■ in the linear case (or when p — > 1/3 in the non-linear case). Unfortunately an 
indeterminate form appears when a 1/2. Is there another way to compute this value? 

On figure[2l obtained by numerical estimations of the integral, the functions p i— > 7p and p i— > 7p 
seem to be convex. In |llj , Volkmer gives a sufficient condition for the convexity of the Lyapunov 
exponent to hold. But this condition is easy to check only in the case of nonnegative matrices, 
and we do not know whether it applies in our cases. 
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